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Effects  Of  Ground  Reflection  On 
Frequency  Modulated  Signals 


1.  INTRODUCTION 

It  is  well-known  that  radar  detection  is  often  adversely  affected 
by  ground  reflections  [Evans.  1966;  Smith  and  Mrstik.  1979] 1,2 . 

Although  several  techniques  have  been  devised  to  alleviate  this  problem 
there  are  still  situations  where  the  problem  can  be  acute.  Particularly 
in  low-angle  radar  tracking,  the  problem  is  sometimes  irremediable 
[Barton.  1974] 3.  The  fact  that  the  ground  surface  is  almost  always 
irregular  introduces  further  complications  and  demands  careful  study 
[Papa  et  al,  1983]  * .  It  must  be  noted  here  that  all  studies  thus  far  on 
the  effects  of  ground  reflections  have  been  restricted  to  conventional 
pulsed  radars.  However,  chirp  radars  use  frequency  modulated  signals  and 
the  effects  of  ground  reflections  on  such  signals  may  be  quite  different 
from  those  on  ordinary  pulsed  signals.  This  report  concentrates  on  the 
effect  of  ground  reflections  on  frequency  modulated  signals. 


2.  DESCRIPTION  OF  THB  PROBLEM 

The  geometry  of  the  problem  is  shown  in  Figure  1.  The  ground  is 
represented  as  a  perfectly  conducting  random  surface  whose  mean 
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FIGURE.  1  SCATTERING  GEOMETRY 


coincides  with  the  xy- plane.  Moreover  the  surface  is  assumed  to  be 
isotropic  and  smoothly  varying  in  slopes  so  that  the  Kirchhoff  method 
may  be  used  for  analysis. 

The  main  quantity  of  interest  to  us  in  this  report  is  the  average 
power  of  the  signal  received  by  the  chirp  radar.  There  are  four 
contributions  to  this  quantity  as  shown  in  Figure  1.  The  target  is 
assumed  to  have  unit  reflectivity.  Both  the  target  and  the  radar  antenna 
are  assumed  to  be  isotropic .  The  ground  reflection  that  we  refer  to  in 
this  report  is  the  specular  reflection  or  the  coherent  reflection. 

Let  fj(t).  f2(t).  f^Ct)  and  f4(t)  represent  the  received  signals 
corresponding  to  the  four  cases  shown  schematically  in  Figure  1 .  The 
angle  between  the  coherent  ray  and  the  z-axis  is  denoted  as  6  while  the 
reflection  coefficient  associated  with  this  ray  is  denoted  as  R.  With 
these  notations  we  have  the  relation 

f2(t)  -  R  f x (t+D)  (1) 

where  D  is  the  delay  caused  by  the  path  difference  between  the  direct 
signal  (Figure  la)  and  the  ground  reflected  signal  (Figure  lb) .  Also  it 
is  clear  that 

f3 (t)  -  f2(t)  (2) 

For  the  same  reason 

f4(t)  -  R2  fx(t+2D)  (3) 

Thus  the  received  signal  f(t)  is  related  to  f^(t)  as 

f(t)  -  fx(t)+  2R  f x (t+D)  +  R2  fx(t+2D)  (4) 

Hereinafter  we  shall  use  capital  letters  to  denote  Fourier  transforms  of 
corresponding  signals.  For  example 
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F^w)  -  J  dt  fx(t)  exp(iwt) 


is  the  Fourier  transform  of  the  signal  f^(t).  Fourier  transformation  of 


Eq.(4)  leads  to 


F(w)  -  [  1  +  2Re  iwD  +  R2e  i2wD  ]  F^w) 


The  signal  transmitted  by  the  chirp  radar  fQ(t)  is  a  linear 
frequency-modulated  pulse  defined  as: 


fQ<t)  -  j  exp[  -  i(wct  +  .5  mt2)] 


for  |t|  *  T 


for  |t|  >  T 


where  wc  is  the  angular  carrier  frequency,  T0  is  half  the  pulse  width, 
and  m  is  the  parameter  determining  the  extent  of  frequency  modulation. 

On  reception  the  signal  passes  through  a  matched  filter  whose 
transfer  function  H(w)  is  given  as 


H(u)  -  exp  [  -  (wc  -  u)2  ] 


This  results  in  an  effective  pulse  compression  [Cook,  1960;  Cohen, 
1987]5,6.  Denoting  the  signal  after  this  filtering  as  g(t)  and  its 
Fourier  transform  as  G(w) , 


G(«)  -  H(u)  F<«) 
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On  substituting  Eq. (6)  into  Eq. (9) 

G(u)  -  Gq(w)  +  Gx(«)  e“iwD+  G2(«)  e"l2wD  (10) 

where 

Gq(w)  -  H(w)  F^(«)  (Ha) 

G^w)  -  2  H(«)  R  F^u)  (lib) 

G2(w)  -  R2  H(w)  F^u)  (11c) 

Inverse  Fourier  transformation  of  Eq.(10)  leads  to 

g(t)  -  gQ(t)  +  g^t+D)  +  g2(t+2D)  (12) 

Notice  that  gQ(t)  is  the  direct  signal  while  gj(t+D)  and  g2(t+2D)  are 
the  ground  reflected  signals  with  delays  D  and  2D  respectively.  Quite 
often  in  practice  these  delays  are  very  small  and  are  much  less  than  the 
pulse  width  of  the  signal.  This  inevitably  results  in  an  interference 
pattern  that  can  be  detrimental  to  radar  detection  and  tracking.  Clearly 
this  problem  is  most  acute  when  g^(t)  or  g2(t)  is  of  the  same  magnitude 
as  gg(t).  For  this  reason,  a  significant  part  of  this  report  will  be 
devoted  to  the  study  of  gj(t)  and  g2(t)  for  various  kinds  of  rough 
surfaces . 

In  Section  3  we  study  the  direct  signal  (Figure  la) .  As  mentioned 
above  there  are  two  ground  reflected  signals  gj^(t)  and  g2(t).  One 
ground  reflection  is  involved  in  g^(t)  whereas  g2(t)  undergoes  two  such 
reflections  (Figure  Id) .  These  reflections  are  indeed  coherent 
reflections  and  they  invariably  depend  on  the  statistical 
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characteristics  of  the  surface  under  consideration.  In  fact  the  coherent 
reflection  coefficient  R  is  given  as  %(v  )  [Beckmann  and  Spizzichino. 
1963] 7  where  x  is  the  characteristic  function  and  vz  =  2  0)  cos0  . 

We  consider  several  statistical  distributions  and  in  each  case  study  the 
characteristics  of  ground  reflected  signals.  The  characteristic 
functions  of  the  four  statistical  distributions  considered  in  this 
report  are  given  in  Appendix  A.  The  corresponding  coherent  reflection 
coefficients  are  given  in  Appendix  B.  Section  4  is  devoted  to  the 
analysis  of  gj(t).  for  each  statistical  distribution,  explicit 
expressions  are  derived  for  g^(t).  In  Section  5  similar  expressions  for 
g2(t)  are  obtained  for  normal  and  exponential  distributions.  The 
results  thus  obtained  are  studied  in  comparison  to  one  another  in 
Section  6 .  The  conclusions  are  presented  in  Section  7 . 


3.  DIRECT  SIGNAL 

Since  we  have  assumed  that  the  target  has  unit  reflectivity,  the 
direct  signal  f1(t)  is  identical  to  the  transmitted  signal  fo(t)  apart 
from  the  time  delay,  which  depends  on  the  range.  For  the  analysis  in 
this  report,  this  time  delay  is  irrelevant.  We  therefore  let  fj(t)  " 


f2(t)  .  Thus 


It  follows  from  Eq.(lla)  that 


g. 


.<*>  '  if  S 


dw  H(w)  F1(«) 


~  iwt 


o 


where 

eg 

Ax  “  /  exP  (“  2^"  <w  ~  2wwc>-  r>) 

“  CO 

-  (2*m  )’5  e  exp  [  (w^-  mt  +  mr)^  j 

On  substituting  Eq. (15)  into  Eq.(14)  and  evaluating  the  integral  we 
obtain 


80<1:)  “  (  2^"  )‘5  exP  [  i(  •5ot2"  wct  "  */4  )  ]  2  T0  sinc  (“To6] 

Thus  the  direct  signal  has  the  anticipated  sharp  peak  at  t  -  0  and  a 
pulse  width  of  -^r- .  Noting  that  the  pulse  width  of  the  transmitted 
signal  is  2Tq  we  see  that  the  ratio  of  pulse  compression  achieved  is 
2mT*  /it . 

O 


(14) 

(15) 

(16) 

4  dB 
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4.  OROOHD  REFLECTED  8X000.  X 

Taking  the  inverse  Fourier  transform  of  Eq.(llb)  ve  get 

*i(t>  -  /  **  2R  H(w>  Fi(w)  e" iwt  (l7) 

—  00 

We  note  that  the  reflection  coefficient  R  for  a  normal  distribution  is 
nondispersive  since  its  statistical  characteristics  are  symmetrical. 
However,  in  applications,  one  notices  that  many  surfaces  do  not  possess 
this  property.  To  study  such  situations  we  consider  as  examples  three 
other  distributions  that  have  asymmetric  characteristics.  The  reflection 
coefficients  R  associated  with  these  statistical  distributions  are  given 
in  Appendix  B.  Below  we  analyze  each  case  for  the  ground  reflected 
signal  Z. 

4.1  Normal  Distribution 

Substituting  Eqs.(8).  (13)  and  (B3)  into  Eq. (17)  we  get 

gi(t)  “  dr  f  **  2  exp  [  -•5°2w2  -  2i"  ( wc~ w )2  ] 

—  <0 


J  dr  exp  ( 
T 


iw  r 
c 


2 

*.5imr  +  iwr 


~  iwt] 


(18a) 


-  J  {  -  5"  "c2  ] 


f  dr  exp  (-  iwcr 
T 


•.5iar2)  A2 


(18b) 
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where  we  have  used  Eq.(B3)  which  gives  R”  exp(-0.  for  the 

reflection  coefficient. 


a  *  2  ( n  e  )°'5  a  cos  0 . 


<e  2 

A2  -  J  dw  exp [  "  (  f“  -  2^-  )  «2  -  i(  t-  r  -  c^/m  )  ] 


“(  */P0  )*5  6Xp[  "  Zt  (  wc/m  +T_  C  ^  ] 


and  where 


K  2  A  i 

pn  *  ^  *  3T 


Substituting  Eq. (19)  back  into  Eq.(18)  and  simplifying,  we  obtain 

*l(t)  -  (  *P0)’'5  “»[  -  ST  "c  -  it  (  t  -  V"  )  2  ]  A3 


where 


A3  " 


/  dr  exp  (  -  a. r  -  b.r) 

-  T  A  L 

o 


im  ,  1 

al  “  2  4pf 


and 


>.  -  »  f  i  +  -J—  ]  - 

1  c  l  2pQm  >  2px 


(19) 


(20) 


(21) 


(22) 


(23a) 


(23b) 
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The  integral  in  Eq.(22)  may  readily  be  expressed  in  terms  of  error 


functions  as  follows 

b2 

Aj  -  0.5  (  ./.j  )  0  5  «cp[  -ji-  ] 

•  {  «f[  Aj  (  V  2^  )  ]  ‘  ,rf[  <24> 

Thus  the  explicit  solution  for  the  ground  reflected  signal  I  for  a 
normal  distribution  of  the  ground  surface  is  contained  in  Eqs.(21).  (23) 
and  (24) . 


4.2  Exponential  Distribution 

From  Eqs . (8) ,  (13).  (B5)  and  (17)  ve  get 

to  iotw 

gi<c>  -  dr  f  *• 2  pis-  “p  I  -  ( v  ")2  ] 

—  CO 

T° 

•  J  dr  exp  (  -  iu  r  -.5imr  -  iwr  -  iwt  ] 

-  T  c 

o 

1  r  i  2  1 
■  x  e*P  I  -  2^  uc  J 

T° 

•  J  dr  exp  (  -  i»cr  -.5imr  ]  (25) 

To 
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vh«r« 
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where 


°1  “  V  "  *  “  °  (32) 

T° 

A,  -  J  dr  exp  (  -  iar(t+a)J 
3  -  T 

o 

exp  [-  2a/9*  ]  erfc  (  17(2*)^  )  (33) 


It  does  not  appear  possible  to  express  in  terms  of  known  functions. 
But  we  note  the  following.  Since  mt~  do  ►  mr  in  our  domain  of 

interest,  namely.  T0  £  t.  if  T  £  T0  it  is  clear  that 

exp  2a/)^  j  erfc  (  L/(2a)^  ]  is  a  slowly  varying  function  of  T.  We 
therefore  approximate  it  by  its  mean  value 

exp  [-  2a/?2  ]  erfc  (  i JTT*)P2  ) 
where 


When  this  approximation  is  used.  Eq.(33)  becomes 


A5  -  exp  [-  2a^2  ]  erfc  (  i JT^Dfi2  )  Ag 


where 

To 

Ag  -  /  dr  exp  (  ~  iar(t+a)] 

-  T 

o 


-  2TQsinc (  mTo(t+a)j 


(34) 


(35) 


(36) 


13 


14 


(40c) 


I3  -  J  d»  exp  [  iw(r-  t)-  ~  (w2-  2uwc)j 


«  mm  »2  2  » 

tjmau  exp  (-  LJkqu  -  a  u  j  erf  (-  law) 


I 2  and  I2  may  be  readily  evaluated  in  terms  of  elementary  functions  as 
follows : 

X1  -  (2*»)‘5  «*p[  "  i*/4  +  i|  (  «c/m  +  r  -  t  “  yira  )2  ]  (4] 


-  b2  r  b2  1 

h  *  “  la*  TTs  exn  4aT~  J 

Z*2  1 

where 


(43a) 


b-  -  -  i(  «  /■  +  r  -  t  -  ^?o  ) 


(43b) 


It  turns  out  that  I3  can  be  expressed  in  terms  of  a  confluent 
hypergeometric  function, 


r2  b 


«  -  JSL—  *  f  2  1.3  3  .  _o_  __2  *J 

1.5  ’l  l  2  '  2  *  2  *  2  ’  a  ’  4a„  J 


2  2 


where 


^(a.bjc.c* ;z,f)  -  ^ 


k-0  t-0 


(*?k+f  (b)k  2k  f* 
(c)k  (c')t  k!  t! 


rfa+kj 
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With  these  results  w«  can  write  g^Ct)  as 

*l(t)  -  i  «p[-  »c2]  (J1+  J2  +  J3]  (47) 

where 


J,  -  /  °dr  exp  (-  i ur  -.5inr2)  I#  ,  t  -  1,2,3 

^  m  C  *  K 

*o 

Jj  is  readily  evaluated  as 

J1  -  (2*b)'5  2T0  exp[  -  iw/4  +  i|  fl2  ]  sine  [  (t+A«  )®tJ 
Oj  “  wc“  mt  -  ajna 

After  some  rearrangements  J2  may  be  written  as 


where 


(48) 


(49) 

(50) 


(51) 


(52) 


and  where 


(53a) 


n„ 


2aom 


+  i« 

c 


(53b) 


a0 


-2 

a 


(53c) 
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On  evaluating  the  integral  in  Eq.(52). 


(54) 


where 


-  exp 


(55a) 


-  [  V 


mb 

2a„ 


Hty 


] }  -  «f{  j. 3  [- 


T  + 
o 


(55b) 


It  is  apparent  that  the  integral  in  J3  is  too  complicated  for  analytic 

evaluation.  But  since  u  /m  -  t  -  ft ro  ►  r  in  -  T  £  r,t  S  T  we  note  that 

c  00 

'Pj  is  a  slowly  varying  function  of  1  in  the  domain  of  the  integrand. 

This  enables  us  to  approximate 


3 

2 


3 

2 


its  mean  value 


Thus 


J 


3 


[3133 

L  2  ’  2  •  2  ’  2 


(56) 
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where 


w  2 

-  exp  [  i  ]  2  (*/a)  ‘5  (1"  i} 

|  erf  [  .5jn  (1+i) (  To+wc/n  )  ] 

-  erf  [  .5jm  (1+i)  (-  To+«e/m  )  ] 

Equation  (47).  together  with  Eqs.(49).  (51).  (54).  (55).  (56).  and  (51) 
provide  the  explicit  solution  for  the  g^(t). 

4.4  Skanr-Noanal  Distribution 

Substituting  Eqs.(8).  (13).  and  (B8)  into  Eq.(17)  we  have 


(57) 


40 

gl(t)  -  “  J  dw  2  exp  [  -  (  «c-  «)2  ] 

•  40 

•  2r exp  [ 1  ( \  l’5  (or  a2)  w  ] 

•  |  exp(  -.Sa^2  u2  )  erfc[  -  ctj«  j 

+  o2  exp(  ~-5o22  w2  )  erfc[  -  o2»  ]  j 


r  °  r  2  a 

J  dr  exp  [  -  i«  r  -.5imr  +  iwr  -  iwt  J 


"  T, 


(58) 
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Let 

8l(t)  "  sl(t)  +  s2(t> 

where 


*1<C)  ‘  2T  /  «“  2  «P  [  -  5-  (  «c-  u)2  ] 

”  to 

•  5T  “»  [ 1  ( ;  )'5  <°r  a2>  “  ] 

•  exp(  -.So^2  w2  )  erfc£  -  a^w  J 


r° 

J  dr 


-  Tr 


f  0 

exp  (  -  i«cr  -.5imr  +  iwr  -  iwt  ] 


s2(t) 


_1 

2k 


I 


dw  2  exp 


2^  [ 1  ( l  )'5  <°r  -2' u  ] 


exp(  “.Sotj2  )  erfc[  a^a  j 


r  °  r  2  •» 

J  dr  exp  [  -  iw  r  ~.5imr  -  Iwr  -  iwt  J 

To 


(59) 


(60a) 


(60b) 
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L«t  us  first  consider  s1(t).  The  RHS  of  Eq. (60a)  may  be  rearranged  as 
follows . 

a\  To  2  \ 

*l(t)  "  UT  J  dr  exp  l  ”  iwc  "•5ittr  J  A8  (61) 

To 

vhere 


A8  -  J  <i“  «*P  [  '  (  V  “)2+  iw(r“  c>] 


•  exp  [  “.Saj2  w2  +i  (  7  )‘5  (“-l-  u  ]  erfc[  “  ~J2~  ] 


*Xp  [  "  2i“  WC2  ]  /  exP  (  “  P]"2"  )  erfc  (  ) 

"  ®0 


and  where 


(62) 


-  .5(  a,2  +i/m)  (63a) 

(63b) 


and 

T1  -  -  lal/j2  (63c) 

The  integral  in  Eq.(62)  may  be  represented  in  terms  of  the 
hypergeometric  function  mentioned  earlier. 
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Substituting  Eq. (64)  into  Eq. (61)  and  rearranging  we  get 


(64) 


sx(C) 


a 


where 


2m 


(65) 


(66) 

(67a) 

(67b) 


(63) 


Note  that  is  the  same  integral  as  in  Eq.(22)  -  Only  the  variables  are 
now  different.  Here  the  result  may  immediately  be  obtained  from  Eq.(24). 
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In  the  integral  in  Eq. (68)  we  nay  employ  the  same  approximation  as 
before,  namely,  identify  the  slowly  varying  part  of  the  integrand  by  its 


mean  value.  To  be  more  specific. 


2  2 
n 


,♦  ra  i.a  i .  „Zi_  _2i  i 

q  1  L  2  ’  2  1  2  '  2  '  Pi  ’  4Pl  J 


_int  ra  i.a  2 . _ V  I 

.  1  u  W1  L  2  •  2  ’  2  ’  2  *  Pl  ’  4px  J 


Thus 


..  ln  Ji.  .  ri  i.a  2  .  _zi - sc  i  A 

1  °  1.5  fl  l  2  '  2  '  2  ’  2  1  p  *  4p,  J  a7 
px  rl  rl 


Equation  (65)  along  with  Eqs.(70)  and  (57)  provides  the  solution  for 
sx(t) . 

We  note  from  Eq. (60)  that  the  integrand  in  s2(t)  is  very  similar  in 
structure  to  that  in  s2(t).  Thus  the  procedure  for  evaluating  s2(t)  is 
identical  to  that  of  sx(t).  We  therefore  omit  the  steps  and  simply 
present  the  results. 


2  bn 


-,(«-•»(  Vi)'-V  “•>[■  5T  “c2‘  £  + 


E{  Ja2  [ 


V^]j  — { 


Act 

2  2 


Ja2  [" 


Ki  ! 


- .« o  w-*  ? *,[  i  ■  i = § .  i 


2  r2 


|  erf  [  .5jm  (1+i)  (  T^/m  ]  ] 


-  erf  [  .5,/m  (1+i)  (-  T+u/m  )  ]  j  (71) 
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where 


a2  -  .25/p2  +  i.5m  (72a) 
t>2  -  .50/p2  +  iu.  (72b) 
P2  "  -5(  «22  +i/»)  (72c) 
72  “  i«2A/2  (72d) 


This  completes  our  derivation  for  g^(t). 

5.  GROUND  REFLECTED  SIGNAL-II 

On  taking  the  inverse  Fourier  transform  of  Eq.(llc) 
i  »*  5  —  lot 

g2(t)  -  ~  J  d«r  H(w)  F1(w)  e  (73) 

""  CO 

In  this  section  we  evaluate  Eq.(73)  for  two  cases  -  normal  distribution 
and  exponential  distribution. 
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Distribution 


Substituting  Eqs.(8).  (13)  and  (B3)  into  Eq.(73)  ve  have 


gjCt)  “  ”2?  _  /  dw  e"  iwt  exp [  -  —  (7-  w)  -  Q2W2  j 


f°  /  2  \ 

•  J  dr  exp  iw  r  -.5imr+iwr  j 


On  comparing  Eq. (74)  with  Eq.(18a)  we  can  immediately  write  down  the 


result  as 


g9(t)  -  .5  g. (t,2o  ) 


5.2  Bspanential  Distribution 


From  Eqs.(8).  (13).  (B5)  and  (73)  we  have 


«o  f  -  iaw  1*5 

g2(t)  -  “2 *„[**{  11757  )  exp[  "  2m”  <V  w>  J 


•  J  dr  exp  iw  r  -.5imr2+iw(r-  t)2/ 


1  r  i  21 

-  exp  r  asr  wc  J 


r°  /  2\ 

J  dr  exp  i»cr  -".5imr  /  B^ 
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where 


m  f  -  low  l2 

1  ■  .  [  **  1  iTJsT  /  •I‘p[  '  2m  <V  2“"c)2+1“<,_  *>] 


-  -§y-  exp  {  .251*  +.25iE2/a  } 


[  (  f  )‘5“  i»^2exp(-  *02)  erfc(i^27«)  ] 


where 


E  -  b  - 


(77) 


(78a) 


(78b) 


Substituting  Eq.(77)  into  Eq. (76)  and  simplifying  we  get 

g2(t)  -  -k-  “P  {  «i»  -  2^  "c  +  -5in3"  /  [  ( ;)'5  V  s2  ]  <79) 


where 


rT°  /  \ 

S.  -  J  dr  exp  \~  lmr(t+2o )y 

^  _  fl* 


(80a) 


-  T 


*2~ 


Jo  /  \ 

-la  J  dr  exp  imr(t+2o)/ 


-  T, 


02  exp(-  a02  )  erfc(l^27a) 


(80b) 


°3  -  °r a 


(80c) 
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Let  us  first  consider  the  integral  in  Eq. (80b) .  Notice  that 
erfc(i02  )  is  a  slowly  varying  function  of  %  and  hence  nay  be 
approximated  by  its  Bean  value  ?2  •XP(~  *?2  )  erfc(ij§2./a)  •  Here 


(81) 


Thus 


*2 - exP(~  «?2  *  •rfc<i/*2 </«)  ^ 


(82) 


^  is  readily  evaluated  as 


i  -  2  T  sincf  ta(t+2ci)T  1 
X  O  v  o  ' 


(83) 


Putting  all  these  together 


g9(t)  -  -§2“  exP  {  -251w  "  2m"  "c  +  *5in3m  ) 


[  (  “)'5~  ia^2  exp("  *^2  *  erfc(1027a)]  3^  (84) 


6.  DISCUSSION 

In  Sections  4  and  5  we  have  calculated  the  ground  reflected 


signal  of  the  chirp  radar  generated  by  four  different  statistical 
distributions.  In  this  section  we  analyze  these  results  and  offer  some 
physical  interpretations.  We  first  look  at  GRS-I. 


C.l  Oramd  Reflected  Signal  -  x 

For  normal  distribution  gj(t)  is  given  by  Eqs.(21)  and  (24). 
Although  it  is  elementary  to  coaputs  this  axprassion  to  obtain  a 
nuaarical  valua.  tha  prasant  fora  doas  not  of far  a  physically 
transparant  pictura.  In  othar  words,  tha  charactaristics  of  tha  rasults 
ara  not  iaaadiataly  avidant  from  tha  axprassions  in  Eqs.(21)  and  (24). 

Ona  way  to  gain  insight  is  to  examine  tha  local  bahavior.  Recall 
that  tha  signal  has  a  peak  at  t  ■  0.  Moreover  our  primary  interest  is 
in  tha  study  of  tha  relative  amplitude  of  the  ground  reflected  signal 
with  respect  to  tha  direct  signal.  We  therefore  look  at  tha  case  whan 
t  ■  0.  Va  further  simplify  our  task  by  restricting  our  attention  to  tha 
region  where  a2a>2  <  1.  Thus  to  a  first  order  in  a2©2 

gx(0)  -  2  (  JS  )•*  xe  [i-.SaV]  (85) 

Note  that  oudc  is  tha  Rayleigh  parameter.  Hence  a2(D2  <  1  refers  to  the 
small  roughness  region.  When  aa>cm‘  0  .  that  is.  when  the  surface  is 
perfectly  smooth.  Eq.(85)  reduces  to 

I  -  2To  (  2m/w)  5  (86) 

Also  from  Eq.(16)  we  see  that 

l8<>(0)|  "  To  (  2a/*)  ’5  (87) 

This  implies  that 

!«!<©> |  -  2  ISqCO) |  (88) 
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froa  Appendix  B  it  is  clear  that  when  cuac  *  0.  R  =  1 ,  hence  it  follows 
that 

gx(t)  “  2g0(t)  (89) 

Thus  it  is  seen  that  Eq. (85)  is  in  agreement  with  Eq.(89).  Also.  Eq.(85) 
suggests  that  |  g^(0)|  should  decrease  with  increasing  OCOC.  This  is  in 
agreement  with  the  physical  fact  that  the  coherent  return  should 
decrease  with  an  increase  in  surface  roughness. 

We  turn  our  attention  now  to  the  exponential  case.  As  before  we 
look  at  the  special  case  when  a  to*  <  1  and  t  “  0 .  It  turns  out  that  to 
a  first  order  in  a20>2 .  |  g1(0)|ejip  has  the  same  expression  as  Eq. (85) . 
This  means  that  to  a  first  order  in  a  .  coherent  returns  from 
normally  distributed  and  exponentially  distributed  surfaces  are 
identical . 

It  is  apparent  that  the  results  obtained  for  surfaces  with  Rayleigh 
and  skew-normal  statistics  have  fairly  complicated  structures.  Even 
from  a  computational  point  of  view  the  resulting  expressions  in  the 
above  cases  are  unmanageable.  It  turns  out  that  for  parameters 
appropriate  for  our  problem  the  series  representation  of  the 
hypergeometric  function  is  not  suitable  for  actual  computations.  Unless 
further  study  of  these  special  functions  are  undertaken  the  value  of 
these  results  may  not  be  fully  appreciated. 

6.2  Oround  Reflected  Signal  -  II 

Consider  the  results  obtained  for  normally  distributed  surfaces . 

When  t  “  0  and  a20>2  <  1  we  have 
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ISjCO)  I  -  To  (  2m/«)  -5  {  1  -  2  aV  } 


(90) 


Notice  that  whan  aa>c  "  0.  Eq.(90)  lands  to 

l«2<0>*  "  !«0<°>l  (91) 

which  agraas  with  tha  physical  fact  that  g2(t)  “  gQ(t)  whan  aa>c  -  0. 
Apart  from  having  all  siailar  charactaristics  as  tha  corresponding 
|gx (0)|  wa  note  from  Eq.  (90)  that  |g2(0)|  decreases  much  faster  with  otoc 
than  |gx (0)|. 

6.3  in— rloal  Results 

Figure  2  is  a  plot  of  the  coaputad  results  for  GRS-I.  Mora 
specifically  it  is  a  plot  of  l  g^(t)/g0(t)(  versus  t.  Tha  parameters 
chosen  are  to  an  extent  arbitrary  but  at  the  saae  time  pertinent  to  our 
problem.  The  solid  lines  and  the  dotted  lines  correspond  respectively  to 
normal  and  exponential  statistics.  The  main  signal  centered  around  t-0 
is  the  compressed  pulse.  The  rest  are  sidelobes.  When  o  -  0.01 
(Figure  2a)  the  normal  and  the  exponential  results  coincide  and 
|  gi  (t)/gn(t)|BIJIT  -  2.  Since  o  -  0.01  corresponds  to  OWc  -  0.015 
both  these  results  are  in  agreement  with  our  predictions.  When 
O  “  0.07  or  ao>c  “  0.104  (Figure  2b)  we  still  see  that  the  normal  and 
exponential  results  coincide.  But  notice  now  that  the  amplitude  of  the 
coherent  signal  is  reduced.  For  further  increase  in  a  we  note  that 
1  gi(t)/g0(t)lnormal  d“r“«"  faster  than  J  gi(t)/g0(t)Jexponential.  The 
point  to  note  is  that  in  some  situations  these  results  can  be  quite 
different  from  one  another. 
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FIGURE  2B 
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FIGURE  2E 


1GHz 


In  Figure  3  wa  have  plotted  GRS-II.  Recalling  the  fact  that  for 
OUDg  "  0.  g2(t)  “  2gj(t)  we  have  appropriately  chosen  the  scale  so  that 
the  results  for  g2(t)  may  be  easily  compared  with  those  of  g^(t).  It  is 
evident  that  g2(t)  has  all  the  characteristics  of  g^t).  But  notice  that 
the  rate  of  decrease  of  |  g2(t)|  with  ouac  is  much  higher  than  that  of 
I  Ct:) |  -  a  fact  clearly  in  agreement  with  our  predictions. 

In  applications,  one  often  models  the  ground  as  a  normally 
distributed  surface.  This  is  primarily  because  such  a  model  facilitates 
easy  analysis.  But  there  is  growing  experimental  evidence  that 
contradicts  this  assumption.  For  example,  terrain  and  a  large  class  of 
sea  surfaces  do  not  have  the  symmetric  roughness  geometry  of  a  normally 
distributed  rough  surface.  To  be  more  precise,  the  valleys  are  often 
more  shallow  and  the  hills  more  steep. 

The  exponential  distribution  is  a  good  candidate  for  an  asymmetric 
rough  surface.  The  results  obtained  in  this  report  bear  out  that  the 
assumption  of  normal  statistics  can  be  meaningful  only  when  cuoc  is 
small.  Otherwise  this  assumption  can  lead  to  significant  errors. 

The  exponential  distribution  is  a  rather  extreme  case  of  asymmetric 
distribution.  A  model  more  suitable  for  application  is  the  Rayleigh 
distribution.  For  terrain,  perhaps  the  most  appropriate  distribution  to 
use  is  the  skew-normal.  Further  study  of  the  results  pertaining  to  these 
models  is  left  for  future  work. 


36 


37 


38 


.2m 


39 


40 


t  (SECONDS) 


41 


FIGURE  3E 


1GHz 
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t  (SECONDS) 


In  this  report,  the  effects  of  ground  reflections  on  frequency 
modulated  signals  were  studied.  On  recognizing  the  detrimental  effect 
on  target  detection  caused  by  interference  due  to  ground  reflections  our 


primary  objective  is  to  compute  the  relative  amplitude  of  the  ground 
reflected  signals  with  respect  to  the  direct  signal.  Four  different 
surfaces  are  considered  for  the  analysis  of  the  ground  reflected 
signals.  Among  them  the  normal  distribution  has  symmetric 
characteristics  while  the  other  three  are  asymmetric .  Due  to  the 
nonlinear  nature  of  the  matched  filter,  the  received  signals  for  various 
types  of  ground  have  quite  different  characteristics.  Notice  that  the 
coherent  reflection  coefficients  of  the  four  surfaces  have  distinct 
frequency  dependences.  Explicit  analytic  expressions  for  the  ground 
reflected  signal  I  are  provided  for  each  of  the  four  types  of  surfaces. 
For  the  sake  of  illustration  we  have  compared  the  ground  reflected 
signals  for  exponentially  distributed  and  normally  distributed  surfaces. 
It  is  observed  that  when  the  Rayleigh  parameter  is  very  small  the  two 
results  are  very  close  to  each  other.  For  larger  Rayleigh  parameters  the 
results  are  quite  different  from  each  other.  This  means  that  the 
universally  accepted  Gaussian  model  may  at  times  lead  to  significant 
errors.  We  conclude  therefore  that  for  the  study  of  scattering  from 
rough  surfaces  a  careful  choice  of  the  statistical  distribution 
appropriate  for  the  problem  is  essential. 
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Appendix  A 

Characteristic  Functions 

In  this  appendix  we  consider  four  statistical  distributions  and 
obtain  their  characteristic  functions.  The  characteristic  function  %(q) 
is  defined  as: 

00 

X<q)  -  J  pco  «xp(iqO  (Al) 

—  00 

where  p(£)  is  the  probability  density  function  of  the  height  of  the 
randoaly  rough  surface.  The  rough  surface  has  zero  mean  and  a  root  mean 
square  height  of  a. 

Al.  NORMAL  DISTRIBUTION 


P(C)  - 


(A2) 


x(q) 


(A3) 
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JU.  KZPOMmiU  OZ8TXXBOTIOH 


~  exp  (  -  ($+a)/<7  )  ;  C  ^  ~ 


;  €  <  ~  o 


x(q)  - 


e~  1qg 

1  ~  iqc 


A3.  RAYLBIQB  DISTRIBUTION! 


[  £  +  ]  e*p|  “  (c/*)2  [  ( 


+  2TJ  / 


for  £  * 


-  J*a 

2c 


otherwise 


where 


c  -  1  -  */4 


{  l  +  i  ^23  etfc  [  - 


i^f]exp( 


-  (,5qo/c)2  ) 
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m.  smr-MOKiAL  DznxzBonoa 


*«>  *  j  77^7  ,xp  (  -  5  («  -  4 

.Trow  *xp  f  -  -5  ((£  •  *)A2)2] :  £  < 

where 

a  “  (2/*)' 5  («t2-  ax) 

<r  -  ^(aj+aj) 

*(q)  -  27"  exp  (iqa) 

|  ax  exp  (-  .5q2«72  )  erfc[  -  ^  ^  ] 

+  tr2  exp  (-  .5q2a2  )  erfc[  -  ^  q*2  ]  J 


(A9) 


(AlOa) 

(AlOb) 


(All) 
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Appendix  B 

Coherent  Reflection  Coefficient 

The  ooherent  reflection  coefficient  R  is  given  by  Beckmann  and 
Spizzichino  [1963] 7  as 

R  *  X(v2)  (Bl) 

where 

vz  -  2  (H0e0)0'5  to  cos  6  (B2) 

Since  the  characteristic  functions  for  various  statistical  distributions 
are  given  in  Appendix  A  the  corresponding  reflection  coefficients  may 
readily  be  obtained  as  follows. 

Bl.  NORMAL  DISTRIBUTION 

R  “  exp  (  -  0.5  a2<02  )  (B3) 

where 

a  -  2  (n0e0)°-5  o  cos  0  (B4) 


B2.  monanu  distribution 


-  iOMO 

s  -  e 

1-  iOCO 


(B5) 
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B3.  RAYLKIGH  DISTRIBUTION 

R  -  exp  (  -  ) 

•  £  1+L/xaw  erfc(iaw)  exp(-  a2u2  ]  j 

where 

o  -  (fie)  -  cos  8 
o  o  c 

B4.  SKEW-NORMAL  DISTRIBUTION 

R  “  2or  exp[  iCW5  (  ar  a2)  «  ] 

*  ^  °y  exp(  “  -5  a2  «2  ]  erfc  [  -  a^w  J 

+  o2  exP(  ”  .5  a2  w2  )  erfc  [  “  ^  a2u  ]  J 

where 

“l  "  2  (Vo),S  °1  cos  9 
*2  “  2  (Vo>'5  °2  cos  # 


(B6) 


(B7 ) 


(B8) 

(B9a) 

(B9b) 


52 


Nomenclature 


» 


%  characteristic  function 
R  reflection  coefficient 

0  delay  in  time  for  the  ground  reflected  signal 
m  frequency  modulation  parameter 
H  transfer  function  of  the  matched  filter 
< a  angular  frequency 
0f>c  angular  carrier  frequency 

2T0  pulse  width  of  the  transmitted  signal 

confluent  hypergeometric  function 
6  angle  of  incidence 
g0(t)  direct  signal 

gj(t)  ground  reflected  signal-I 

g2(t)  ground  reflected  signal -II 
o  rms  height  of  the  rough  surface 
a  2  ( )0’5  a  cos  6 

a  0.5  a  (l  -O^Sjc)-0*5 

a±  2(noeo)°‘50i  cos  0  :  i  -  l,  2 
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ufisds  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other 
ESD  elements  to  perform  effective  acquisition  of  C3I  systems.  In  addition , 
Rome  Laboratory?*  technology  supports  other  AFSC  Product  Divisions ,  the 
Air  Force  user  community,  and  other  DOD  and  non-DOD  agencies.  Rome 
LdB&rai&ry  maintains  technical  competence  and  research  programs  in  areas 
Including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
Li  ul  software  productbillty,  wide  area  surveillance/sensors,  signal  proces¬ 
sor®,  solid  state  sciences,  photonics,  electromagnetic  technology,  super - 
conductivity,  and  electronic  reliability/maintainability  and  testability. 


